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Abstract. We study the problem of generating monomials of a polynomial in 
the context of enumeration complexity. In this setting, the complexity measure 
is the delay between two solutions and the total time. We present two new 
algorithms for restricted classes of polynomials, which have a good delay and 
the same global running time as the classical ones. Moreover they are simple to 
describe, use little evaluation points and one of them is parallelizable. 
We introduce three new complexity classes, TotalPP, IncPP and DelayPP, 
which are probabilistic counterparts of the most common classes for enumeration 
problems, hoping that randomization will be a tool as strong for enumeration 
as it is for decision. Our interpolation algorithms proves that a lot of interesting 
problems are in these classes like the enumeration of the spanning hypertrees 
of a 3-uniform hypergraph. 

Finally we give a method to interpolate a degree 2 polynomials with an accept- 
able (incremental) delay. We also prove that finding a specified monomial in a 
degree 2 polynomial is hard unless RP = NP. It suggests that there is no algo- 
rithm with a delay as good (polynomial) as the one we achieve for multilinear 
polynomials. 



1 Introduction 

Enumeration, the task of generating all solutions of a given problem, is an interesting 
generalization of decision and counting. Since a problem typically has an exponential 
number of solutions, the way we study enumeration complexity is quite different from 
decision. In particular, the delay between two solutions and the time taken by an 
algorithm relative to the number of solutions seem to be the most considered complexity 
measures. In this paper, we revisit the famous problem of polynomial interpolation, that 
is to say finding the monomials of a polynomial from its values, with these measures 
in mind. 

It has long been known that a finite number of evaluation points is enough to 
interpolate a polynomial and efficient procedures (both deterministic and probabilistic) 
have been studied by several authors [11213] . The complexity depends mostly on the 
number of monomials of the polynomial and on an a priori bound on this number 
which may be exponential in the number of variables. The deterministic methods rely 
on prime numbers as evaluation points, with the drawback that they are very large. 
The probabilistic methods crucially use the Schwarz-Zippel lemma, which is also a tool 
in this article, and efficient solving of particular linear systems. 

As a consequence of a result about random efficient identity testing [4] , Klivans and 
Spielman give an interpolation algorithm, which happens to have an incremental delay. 



In this vein, the present paper studies the problem of generating the monomials of a 
polynomial with the best possible delay. In particular we consider natural classes of 
polynomials such as multilinear polynomials, for which we prove that interpolation can 
be done efficiently. Similar restrictions have been studied in other works about identity 
testing (the decision version of interpolation) for a quantum model [5] or for depth 3 
circuits which thus define almost linear polynomials [BJ. Moreover, a lot of interesting 
polynomials are multilinear like the Determinant, the Pfaffian, the Permanent, the 
elementary symmetric polynomials or anything which may be defined by a syntactically 
multilinear arithmetic circuit. 

In Sec. 0] we present an algorithm which works for polynomials such that no two of 
their monomials use the same set of variables. It is structured as in [4] but is simpler and 
has better delay, though polynomially related. In Sec. [5] we propose a second algorithm 
which works for multilinear polynomials; it has a delay polynomal in the numberof 
variables, which makes it exponentially better than the previous one and is also easily 
parallelizable. In addition both algorithms enjoy a global complexity as good as the 
algorithms of the literature, are deterministic for monotone polynomials and use only 
small evaluation points making them suitable to work over finite fields. 

We describe in Sec. |6j three complexity classes for enumeration, namely TotalP, 
IncP, DelayP which are now commonly used |7|8|9|10j to formalize what is an effi- 
ciently enumerable problem. We introduce probabilistic variants of these classes, which 
happen to characterize the enumeration complexity of the different interpolation algo- 
rithms. Their use on polynomials computable in polynomial time enable us to prove 
that well-known problems are in these classes. Those problems already have better 
enumeration algorithms except the last, enumeration of the spanning hypertrees of 
a 3-uniform hypergraph, for which our method gives the first efficient enumeration 
algorithm. 

In the last section we discuss how to combine the two algorithms we have presented 
to interpolate degree 2 polynomials with incremental delay. We also prove that the 
problem of finding a specified monomial in a degree 2 polynomial is hard by encoding a 
restricted version of the hamiltonian path problem in a polynomial given by the Matrix- 
Tree theorem (see Thus there is no polynomial delay interpolation algorithm for 
degree 2 polynomials similar to the one for degree 1 because it would solve the later 
problem and would imply RP = NP. Finally we compare our two algorithms with 
several classical ones and show that they are good with regard to parameters like 
number of calls to the black box or size of the evaluation points. 

Please note that most proofs are given in the appendix. 

2 Enumeration Problems 

In this section, we recall basic definitions about enumeration problems and complexity 
measures and we introduce the central problem of this article. 

The computation model is a RAM machine as defined in [10] which has, in addition 
to the classical definition, an instruction Write(A) which outputs the content of the 
register A. The result of a computation of a RAM machine is the sequence of integers 
which were in A when the instructions Write(A) were executed. For simplicity we 
consider that these integers encode words, and that the input of the machine is also a 
word represented by suitable integers in the input registers. Let M be such a machine 



and x a word, we write M(x) the result of the computation of M on x. The order 
in which the outputs are given does not matter, therefore M(x) will denote the set of 
outputs as well as the sequence. We choose a RAM machine instead of a Turing machine 
since it may be useful to deal with an exponential amount of memory in polynomial 
time, see for instance the enumeration of the maximal independent sets of a graph [7] . 

Definition 1 (Enumeration Problem). Let A be a polynomially balanced binary 
predicate, i.e. A{x,y) =>• \y\ < Q(\x\), for a certain polynomial Q. We write A(x) for 
the set of y such that A(x,y). We say that a RAM machine M solves the enumeration 
problem associated to A, Enum-A for short, if M(x) = A(x) and there is no repetition 
of solutions in the computation. 

Let T(x, i) be the time taken by a machine M to return i outputs from the instance 
x. As for decision problems, we are interested by the total time taken by M, namely 
T(x, \M{x)\). We are also interested by the delay between two solutions, that is to say 
T(x,i + 1) — T{x,i). M has an incremental delay when it is polynomial in \x\ and i, 
and M has a polynomial delay when it is polynomial in \x\ only. 

A probabilistic RAM machine has a special instruction rand which writes in a 
specific register the integer or 1 with equal probability. All outcomes of the instruction 
rand during a run of a RAM machine are independent. 

Definition 2 (Probabilistic enumeration). We say that the probabilistic RAM ma- 
chine M solves ENUM-^4 with probability p if P[A(x) = M{xj\ > p and there is no 
repetition of solutions in the computation. 

We adapt the model to the case of a computation with an oracle, by a special 
instruction which calls the oracle on a word contained in a specific register and then 
writes the answer in another register in unit time. 

In this article we interpret the famous problem of interpolating a polynomial given 
by a black box as a enumeration problem. It means that we try to find all the mono- 
mials of a polynomial given by the number of its variables and an oracle which allows 
to evaluate the polynomial on any point in unit time. This problem is denoted by 
ENUM-PoZy but will be solved in this article only on restricted classes of polynomials. 

3 Finding one Monomial 

In this section we introduce all the basic tools we need to build interpolation algorithms. 
One consider polynomials with n variables and rational coefficients. A sequence of n 
positive integers e = (ex, ... , e n ) characterizes the monomial X e = X^ 1 A| 2 . . . ■ 
We call t the number of monomials of a polynomial P written P(X) = ^jX ej ■ 

}<3<t 

The degree of a monomial is the maximum of the degrees of its variables and the 
total degree is the sum of the degrees of its variables. Let d (respectively D) denote 
the degree (respectively the total degree) of the polynomial we consider, that is to say 
the maximum of its monomial's degree (respectively total degree). In Sec. [5] we assume 
that the polynomial is multilinear i.e. d = 1 and D is thus bounded by n. 

We assume that the maximum of the bitsize of the coefficients appearing in a 
polynomial is 0(n) to simplify the statement of some results, in the examples of Sec. [6] 
it is even 0(1). When analyzing the delay of an algorithm solving ENVM-Poly we are 



interested in both the number of calls to the black box and the time spent between 
two generated monomials. We are also interested in the size of the integers used in the 
calls to the oracle, since in real cases the complexity of the evaluation depends on it. 

The support of a monomial is the set of indices of variables which appears in the 
monomial. Let L be a set of indices of variables, for instance a support, then //, is the 



Prom now on, we denote Jl (P) by Pj, . It is the polynomial obtained by substituting 
to every variable of index not in L, that is to say all the monomials of P which have 
their support in L. We call X L the multilinear term of support L, which is the product 
of all Xi with i in L. 

Lemma 1. Let P be a polynomial without constant term and whose monomials have 
distinct supports and L a minimal set (for inclusion) of variables such that Pl is not 
identically zero. Then there is an integer X such that Pl = XX L . 

The hrst problem we want to solve is to decide if a polynomial given by a black box 
is the zero polynomial, a problem called Polynomial Identity Testing. We are especially 
interested in the corresponding search problem, i.e. giving explicitly one term and its 
coefficient. Indeed, we show in Sec. [4] how to turn any algorithm solving this problem 
into an incremental interpolation algorithm. 

It is easy to see [2] that a polynomial with t monomials has to be evaluated in t 
points to be sure that it is zero. If we do not have any a priori bound on t, then we 
must evaluate the polynomial on at least (d + 1)™ n-tuples of integers to determine 
it. As we are not satisfied with this exponential complexity, we introduce probabilistic 
algorithms, which nonetheless have a good and manageable bound on the error. 

Lemma 2 (Schwarz-Zippel [12j ) . Let P be a non zero polynomial with n variables 
of total degree D, if x%, . . . , x n are randomly chosen in a set of integers S of size — 
then the probability that P(xi, . . . ,x n ) = is bounded by e. 

A classical probabilistic algorithm to decide if a polynomial P is identically zero 
can be derived from this lemma. It picks x\, . . . ,x n randomly in [-^0 and calls the 
oracle to compute P(xi, . . . , x n ). If the result is zero, the algorithm decides that the 
polynomial is zero otherwise it decides that it is non zero. Remark that the algorithm 
never gives a false answer when the polynomial is zero. The probability of error when 
the polynomial is non zero is bounded by e thanks to Lemma [2] Polynomial Identity 
Testing is thus in the class recognizable by a polynomial time algorithmRP. 

This procedure makes exactly one call to the black box on points of size log(-^-)- 
The error rate may then be made exponentially smaller by increasing the size of the 
points. There is an other way to achieve the same reduction of error. Repeat the 
previous algorithm k times for e = \, that is to say the points are randomly chosen 
in [2D]. If all runs return zero, then the algorithm decides that the polynomial is zero 
else it decides it is non zero. The probability of error of this algorithm is bounded by 
2 _fc , thus to achieve an error bound of e we have to set k — log(i). We denote by 
not_zero{P 1 e) the latter procedure, which is given as inputs a black box polynomial P 
and the maximum probability of failure e. It uses slightly more random bits but it only 
involves numbers less than 2D. 

1 We write [x] for the set of integers between 1 and \x] . 



homomorphism of Q[Xi, . . . , X n ] defined 




Up to Sec. [5j all polynomials have monomials with distinct supports and no con- 
stant term. This class of polynomials contains the multilinear polynomials but is much 
bigger. Moreover being without constant term is not restrictive since we can always 
replace a polynomial by the same polynomial minus its constant term that we compute 
beforehand by a single oracle call to P(0, . . . , 0). 

We now give an algorithm which finds a monomial of a polynomial P, in randomized 
polynomial time thanks to the previous lemmas. In this algorithm, L is a set of indices 
of variables and i an integer used to denote the index of the current variable. 



Algorithm 1: find_monomial 
Data: A polynomial P with n variables and the error bound e 
Result: A monomial of P 
begin 

L<— {l,...,n} 
if not_zeroff ,— rr) then 
for i = 1 to n do 

if not-zero (P L \ {i} , 7 ^ T ) then 
[ L L<—L\{z} 

return The monomial of support L 
else 

L return "Zero" 



Once a set L is found such that Pl is a monomial XX e , we must compute A and 
e. The evaluation of Pl on (1, . . . , 1) returns A. For each i 6 L the evaluation of Pl 
on Xi = 2 and for j ^ i, Xj = 1 returns X2 ei . From these n calls to the black box, we 
compute e in linear time and thus output XX e . 

We analyze this algorithm, assuming first that the procedure not_zero never makes 
a mistake. We also assume that P is not zero, which means that the algorithm has 
not answered "Zero" . In this case at the end of the algorithm, Pl is not zero. In fact 
we remove an element from L only if this condition is respected. As removing another 
element from L would make Pl zero by construction, the set L is minimal for the 
property of Pl being non zero. Then by Lemma [T] we know that Pl is a monomial of 
P, which allows us to output it as previously explained. 

Errors only appear in the procedure not_zero with probability ^rj- Since we use 
this procedure n+l times we can bound the total probability of error by e. The total 
complexity of this algorithm is 0(nlog(^-)) since each of the n calls to the procedure 
not-ze.ro makes 0(log(^)) calls to the oracle in time O(l). We summarize the properties 
of this algorithm in the next proposition. 



Proposition 1. Given a polynomial P as a black box, whose monomials have dis- 
tinct supports, Algorithm^ finds, with probability 1 — e, a monomial of P by making 
0(n log(^)) calls to the black box on entries of size log(2P). 



4 An Incremental Algorithm for Polynomials with Distinct 
Supports 

We build an algorithm which enumerates the monomials of a polynomial incrementally 
by using the procedure findjmonomial defined in Proposition [TJ Recall that incre- 
mentally means that the delay between two consecutive monomials is bounded by a 
polynomial in the number of already found monomials. 

We need a procedure subtract(P, Q) which acts as a black box for the polynomial 
P — Q when P is given as a black box and Q as an explicit set of monomials with their 
coefficients. Let D be the total degree of Q, C a bound on the size of its coefficients 
and i be the number of its monomials. One evaluates the polynomial subtract(P, Q) 
on points of size m as follows: 

1. compute the value of each monomial of Q in time 0(D max(C, m)) 

2. add the values of the i monomials in time 0(iD max(C, m)) 

3. call the black box to compute P on the same points and return this value minus 
the one we have computed for Q 



Algorithm 2: Incremental computation of the monomials of P 
Data: A polynomial P with n variables and the error bound e 
Result: The set of monomials of P 
begin 

while not_zero(subtract('P,Q y ), ^ft) do 

M < — find_monomial(subtract(P,Q), ^ft) 
Write(M) 
|_ Q < — Q + M 



Theorem 1. Let P be a polynomial whose monomials have distinct supports with n 
variables, t monomials and total degree D. Algorithm^ computes the set of monomials 
of P with probability 1 — e. The delay between the i th and i + 1 th outputted monomials 
is bounded by 0(iDn 2 (n + log(i))) in time and 0{n{n + log(-))) calls to the oracle. 
The algorithm performs 0(tn(n + log(i))) calls to the oracle on points of size \og(2D). 



5 A Polynomial Delay Algorithm for Multilinear Polynomials 

In this section we introduce an algorithm which enumerates the monomials of a multi- 
linear polynomial with a polynomial delay. This algorithm has the interesting property 
of being easily parallelizable, which is obviously not the case of the incremental one. 

Let P be a multilinear polynomial with n variables of total degree D. Let L\ and 
1/2 be two disjoint sets of indices of variables and I the cardinal of L-2.. We can write 
P LlUL . 2 = X L ' 2 Pi(X) + P 2 {X), where X L<2 does not divide P 2 (X). We want to decide 
if there is a monomial of P, whose support contains L 2 and is contained in L\ U L 2l 



which is equivalent to deciding wether P\{X) is not the zero polynomial. To do this, 
we define a univariate polynomial H(Y) from Pl 1 ul 2 '- 

1. substitute a randomly chosen value Xi in [2D] to Xj for all i G L\ 

2. substitute the variable Y to each Xi with i e I2 

The polynomial H{Y) can be written Y Pi (a;) + Pt,(x, Y). If Pi is a non zero polyno- 
mial then Pi (a:) is a non zero constant with probability at least \ because of Lemma 
O Moreover P 2 (x,Y) is a polynomial of degree strictly less than I. Hence, to decide if 
the polynomial Pi is not zero, we have to decide if H(Y) is of degree I. 

To this aim we do a univariate interpolation of H(Y): for this we need to make I 
oracle calls on values from 1 to I. The time needed to do an interpolation thanks to 
these values, with s a bound on the size of H{i) for 1 < i < I, is 0(/ 2 log(s)). We 
improve the probability of error of the described procedure from \ to e by repeating it 
log(-) times and name it not_zero-improved(Li,L 2 ,P 7 e). 

We now describe a binary tree which contains informations about the monomials 
of P. The set of node of this tree is the pairs of list (Li,L 2 ) such that there exists 
a monomial of support L in P with L 2 C L C L\ U L 2 . Consider a node labeled by 
(£-1,1/2), we note i the smallest element of L\, it has for left child (Li \ {i},L 2 ) and 
for right child {L\ \ {i}, L 2 U {«}) if they exist. The root of this tree is ([n], 0) and the 
leaves are of the form (0,L 2 )- There is a bijection between the leaves of this tree and 
the monomials of P: a leaf (0, L 2 ) represents the monomial of support L 2 . 

To enumerate the monomials of P, Algorithm [3] does a depth first search in this 
tree using not_zero_improved and when it visits a leaf, it outputs the corresponding 
monomial thanks to the procedure coefficient (P, L) that we now describe. We have L 
of cardinality I the support of a term and we want to find its coefficient. Consider H(Y) 
built from L\ — and L 2 = L, the coefficient of Y l in this polynomial is the coefficient 
of the monomial of support L. We interpolate H(Y) with I calls to the oracle as before 
and return this coefficient. 



Algorithm 3: A depth first search of the support of monomials of P, recursively 

written 

Data: A multilinear polynomial P with n variables and the error bound e 

Result: All monomials of P 

begin 

Monomial(Li , L2, i) = 
if i = n + 1 then 

I Write(coefficient(P,L 2 )) 
else 

if not_zero_improved(Li \ {i}, L2, P, 2^7) then 
J_ Monomial(Li \ {«}, L 2 ,i + 1) 

if not_zero_improved(Li \ {i}, L2 U {i}, P, ofc) then 
|_ Monomial(Li \ {«}, L 2 U {i}, i + 1) 

in Monomial([n], 0, 1) 



Theorem 2. Let P be a multilinear polynomial with n variables, t monomials and 
total degree D. Algorithm^ computes the set of monomials of P with probability 1 — 



e. The delay between the i and i + 1 outputted monomials is bounded in time by 
0(D 2 n 2 log(n)(n+log(i))) and by 0(nD(n+log(-))) oracle calls. The whole algorithm 
performs 0(tnD(n + log(^))) calls to the oracle on points of size 0(log(D)). 

There is a possible trade-off in the way not-Zero -improved and coefficient are im- 
plemented: if one knows a bound on the size of the coefficients of the polynomial and 
use exponentially bigger evaluations points then one needs only one oracle call. The 
number of calls in the algorithm is then less than tn which is close to the optimal 2t. 

Remark that when a polynomial is monotone (coefficients all positive or all negative) 
and is evaluated on positive points, the result is zero if and only if it is the zero 
polynomial. Algorithms [5] and [3] may then be modified to work deterministically for 
monotone polynomials with an even better complexity. 

Moreover both algorithms work for polynomials over Q but we can extend them to 
work over finite fields. Since they only use evaluation points less than 2D, polynomials 
over any field of size more than 2D can be interpolated with very few modifications, 
which is good in comparison with other classical algorithms. 

6 Complexity Classes for Enumeration 

In this part the results about interpolation in the black box formalism are transposed 
into more classical complexity results. We are interested in enumeration problems de- 
fined by predicates A(x, y) such that there is for each x a polynomial P x whose monomi- 
als are in bijection with A(x). If P x is efficiently computable, an interpolation algorithm 
gives an effective way of enumerating its monomials and thus to solve Enum-A 

Example 1. We associate to each graph G the determinant of its adjacency matrix. 
The monomials of this multilinear polynomial are in bijection with the cycle covers 
of G. Hence the problem of enumerating the monomials of det(M) is equivalent to 
enumerating the cycle covers of G. 

The specialization of different interpolation algorithms to efficiently computable 
polynomials naturally correspond to three "classical" complexity classes for enumer- 
ation and their probabilistic counterparts. We present several problems related to a 
polynomial as in Example [T] to illustrate how easily the interpolation methods de- 
scribed in this article produce enumeration algorithms for combinatorial problems. 
Although the first two examples already had efficient enumeration algorithms, the last 
did not, which shows that interpolation methods can bring new results in enumeration 
complexity. 

In all the following definitions, we assume that the predicate which defines the 
enumeration problem is decidable in polynomial time, that is to say the corresponding 
decision problem is in P. 

Definition 3. A problem Enum-j4 is decidable in polynomial total time TotalP (resp. 
probabilistic polynomial total time TotalPPj if there is a polynomial Q(x,y) and a 
machine M which solves Enum-A (resp. with probability greater than ^) and satisfies 
for all x, T(x, \M(x)\) < Q(\x\, \M(x)\). 

TotalPP is very similar to the class BPP for decision problems. By repeating 
a polynomial number of times an algorithm working in total polynomial time and 



returning the set of solutions we find in the majority of runs, we decrease exponentially 
the probability of error. The choice of | is hence arbitrary, everything greater than 
| would do. This property holds for the other probabilistic classes we are going to 
introduce, but unlike TotalPP the predicate which defines the enumeration problem 
needs then to be decidable in polynomial time 

Early termination versions of Zippel's algorithm 2 3 solve enumPoly in a time 
polynomial in the number of monomials. If we now use this algorithm on the De- 
terminant which is computable in polynomial time, we enumerate its monomials in 
probabilistic polynomial total time. Thanks to Example [TJ the enumeration of the 
cycle covers of a graph is in TotalPP. 

Definition 4. A problem Enum-A is decidable in incremental polynomial time IncP 
(resp. probabilistic polynomial total time IncPPJ if there is a polynomial Q(x,y) and 
a machine M which solves Enum-A (resp. with probability greater than |j and satisfies 
for all x, T(x, i + 1) — T(x, i) < Q(\x\, i). 

The classes IncP and IncPP can be related to the following search problem, 
parametrized by a polynomially balanced predicate A. 

AnotherSolution^ 

Input: An instance x of A and a subset S of A(x) 

Sortie: An element of A(x) \ S or a special value if A(x) — S 

It has been proved [8] that AnotherSolution^ g FP if and only if A e IncP. 
We adapt this result to the class IncPP. If A is a polynomial predicate, the search 
problem is to return for all x an element of A(x) or a special value if A(x) is empty. A 
search problem has a solution in probabilistic polynomial time if there is a polynomial 
time algorithm which solves the search problem with probability | . 

Proposition 2. AnotherSolution^ has a solution in probabilistic polynomial time 
if and only if A £ IncPP . 

Since Zippel's algorithm finds all monomials in its last step, it seems hard to turn it 
into an incremental algorithm. On the other hand Algorithm [5] whose design has been 
inspired by Proposition [2] does the interpolation with incremental delay. 

Example 2. To each graph we associate a polynomial PerfMatch, whose monomials 
represent the perfect matchings of this graph. For graphs with a "Pfaffian" orientation, 
such as the planar graphs, this polynomial is related to a Pfaffian and is then efficiently 
computable. Moreover all the coefficients of this graph are positive, therefore we can 
use Algorithm [2] to interpolate it deterministically with incremental delay. We have 
then proved that the enumeration of perfect matching is in IncP. 

Definition 5. A problem Enum-A is decidable in polynomial delay DelayP, (resp. 
probabilistic polynomial delay time DelayPP,) if there is a polynomial Q(x,y) and a 
machine M which solves Enum-^4 (resp. with probability greater than ^) and satisfies 
for all x, T(x, i + 1) — T(x, i) < Q(\x\). 

Example 3 (Spanning Hypertrees) . The notion of a spanning tree of a graph has several 
interesting generalizations to the case of hypergraphs. Nevertheless deciding if there is 



a spanning hypertree is polynomially computable only for the notion of Berge acyclic- 
ity and 3-uniform hypergraphs 13 thanks to an adaptation of the Lovasz matching 
algorithm in linear polymatroids |14j . 

A polynomial Z is defined for each 3-uniform hypergraph [15 with coefficients — 1 
or 1, whose monomials are in bijection with the spanning hypertrees of the hypergraph. 
A new Matrix- Tree theorem [TS] shows that Z is the Pfaffian of a matrix, whose co- 
efficients are linear polynomials depending on the hypergraph. Thus Z is efficiently 
computable by first evaluating a few linear polynomials and then a Pfaffian. This has 
been used to give a simple RP algorithm |16j to decide the existence of a spanning 
hypertree in a 3-uniform hypergraph. 

If we use Algorithm [3] we can enumerate the monomials of Z with probabilistic 
polynomial delay. The delay is good since the total degree of the monomials is low 
and the size of the coefficients is 1, which helps in the interpolation of the univariate 
polynomials. As a conclusion, the problem of enumerating the spanning hypertrees of 
a 3-uniform hypergraph is in DelayPP. 

7 Degree 2 Polynomials 

7.1 An Incremental Algorithm for Degree 2 Polynomials 

We now give an incremental algorithm for the case of polynomials of degree d = 2. 
It is enough to describe a procedure which finds a monomial of a polynomial P, then 
Algorithm [2] turns it into an incremental algorithm. 

First remark that we may use Algorithm [I] on a polynomial P of arbitrary degree 
to find a minimal support L in P. Since it is minimal, all monomials of Pl have L as 
support and Pl(X) = X L Q(X) with Q a multilinear polynomial. Therefore if we find 
a monomial of Q(X) and multiply it by X L , we have a monomial of P. 

We may simulate an oracle call to Q(X) by a call to the oracle giving Pj, and 
a division by the value of X L as long as no Xi is chosen to be 0. Remark that the 
procedure not_zero-improved(L' , L\L' , Q, e) calls the black box only on strictly positive 
values since L = L' U {L \ L') . It allows us to decide if Q has a monomial whose support 
contains L' . In Algorithm 0] we find a L' such that it is contained in the support of a 
monomial and is maximal for this property. Since Q(X) is multilinear there is only one 
monomial of support L' and we find its coefficient by the procedure coefficient (Q ,L' ). 

Algorithm 4: Finding a monomial of a degree two polynomial 

Data: A polynomial P L = X L Q(X) of de gree 2 with n variables, an error 

bound e 
Result: A monomial of Q 
begin 

L' i — 

for i = 1 to n do 

if not_zero_improved(0, L' U {i}, Q,-) then 

L /.' • /.' M 

return coefficient (Q, L' ) 

Thanks to Algorithm 2] we have a monomial of Q and if we multiply it by X L it 
is a monomial of P. We then use it to implement find-monomial in Algorithm [1] and 
obtain an incremental interpolation algorithm for degree 2 polynomials. 



7.2 Limit to the Polynomial Delay Approach 

Here we study the problem of deciding if a monomial has coefficient zero in a polyno- 
mial. In the case of a multilinear polynomial the procedure notjzeroJimproved solves 
the problem in polynomial time but for degree 2 polynomials we prove it is unpossible 
unless RP = NP. Therefore there is no generalization of Algorithm [3] to higher degree 
polynomials, although a polynomial delay algorithm may exist. 

Proposition 3. Assume there is an algorithm which, given a polynomial of degree 2 
and a monomial, can decide in probabilistic polynomial time if the monomial appears 
in the polynomial then RP = NP. 

Proof. Let G be a directed graphs on n vertices, the Laplace matrix L(G) is defined 
by L(G)i,j = -X^j when e E(G), L(G) iti = V" X itj and otherwise. The 

(ij)es(G) 

Matrix- Tree theorem is the following equality where T s is the set of spanning trees of G 
whose all edges are oriented away from the vertex s and L(G) s j is the minor of L(G) 
where row s and column t have been deleted: 

det(L(G) M )(-l) s+t = E II X U 

Te% (i,j)£T 

We substitute to Xi.j the product of variables YiZj in the polynomial det(L(G) Sj t) 
which makes it a polynomial in In variables still computable in polynomial time. Every 
monomial represents a spanning tree whose maximum outdegree is the degree of the 
polynomial. We assume that every vertex of G has indegree and outdegree less or equal 
to 2 therefore det(L(G) S! i) is of degree 2. 

Remark now that a spanning tree, all of whose vertices have outdegree and indegree 
less or equal to 1 is an Hamiltonian path. Therefore G has an Hamiltonian path begin- 
ning by s and finishing by a vertex v if and only if det(L(G) Si i) contains the monomial 

Y S Z V n#«,u 

There is only a polynomial number of pairs (s, v), thus if we assume there is a prob- 
abilistic polynomial time algorithm to test if a monomial is in a degree 2 polynomial, 
we can decide in probabilistic polynomial time if G of outdegree and indegree at most 
2 has an Hamiltonian path. Since this problem is NP complete [17] we have RP = NP. 

8 Conclusion 

Let us compare our method to three classical interpolation algorithms, which unlike 
our method can interpolate polynomials of any degree. Once restricted to multilinear 
polynomials, Algorithm |3] is really efficient compared to the algorithm of Klivans and 
Spielman (KS), which is the only known method with a bound on the delay. Note also 
that Algorithm [2 which is not presented in the next table, needs n 2 calls to the black 
box to guess one monomial, whereas KS needs (nD) e calls and then the same method 
is used to recover the whole polynomial from this procedure. 

In the table T is a bound on t the number of monomials that Ben-Or Tiwari and 
Zippel algorithms need to do the interpolation. In the row labeled Enumeration is 
written the kind of enumeration algorithm the interpolation method gives when the 
polynomial is polynomially computable. 





Ben-Or Tiwari [1 


Zippel 2 


KS [4] 


Algorithm [3] 


Algorithm type 


Deterministic 


Probabilistic 


Probabilistic 


Probabilistic 


Number of calls 


2T 


tnD 


t(nDf 


tnD{n + log(i)) 


Total time 


Quadratic in T 


Quadratic in t 


Quadratic in t 


Linear in t 


Enumeration 


Exponential 


TotalPP 


IncPP 


DelayPP 


Size of points 


Tlog(n) 


log(^) 


log(2£) 


Iog(£>) 
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Appendix 



Here we give most of the proofs which are omitted in the article and the alternate 
implementation of not_zero-improved with only one oracle call. 

Proof of Theorem [T] : 
Correction : 

We analyze this algorithm under the assumption that the procedures not-zero and 
find-monomial do not make mistakes. 

We have the following invariant of the while loop : Q is made from a subset of 
the monomials of P. It is true at the beginning because Q is zero. Assume that 
Q satisfies this property at a certain point of the while loop, since we know that 
not-zero (subtract (P ,Q)) , P — Q is non zero and is then a non empty subset of the 
monomials of P. The outcome of find-monomial (subtract (P,Q)) is thus a monomial of 
P which is not in Q, therefore Q plus this monomial still satisfies the invariant. Remark 
that we have also proved that the number of monomials of Q is increasing by one at 
each step of the while loop. The algorithm must then terminate after t steps and when 
it does not-zero (subtract (P ,Q)) gives a negative answer meaning that Q = P. 
Probability of error : 

The probability of failure is bounded by the sum of the probabilities of error coming 
from not-zero and find-monomial. We both call these procedures t times with an error 
bounded by ^tt- Since 2t < 2" +1 , the total probability of error is bounded by e. 
Complexity : 

The procedure not-zero is called t times and uses the oracle n + log(i) times, whereas 
find-monomial is called t times but uses n(n + log(-)) oracle calls, which adds up to 
t(n + l)(n + log(i)) calls to the oracle. In both cases the evaluation points are of size 
0(D). 

The delay between two solutions is bounded by the evaluation of find-monomial, 
which is dominated by the execution of subtract(P, Q) at each oracle call on points of 
size D. The algorithm calls subtract(P, Q) n(n + log(i)) times and each of these calls 
needs 0(iD max(C, £>)), which gives a delay of 0(iD max(C, D)n(n + log(f))). □ 

Alternate method to implement not-zero -improved : 

We want to decide if P\(X) is the zero polynomial in Pl 1 ul 2 — X L2 Pi(X) + P 2 (X). 
We let a be the integer 2 2 ("+ c '+ Dlo s( ; T i )) and I the cardinal of L 2 . We do a call to the 
oracle on the values (xi) ie ^ n y. 

Xi is randomly chosen in [— ] if i G L\ 
Xi — a if i G L 2 

Xi = otherwise 

The value of a variable which is not in L 2 is bounded by — , therefore a monomial 
of P 2 (which contains at most I — 1 variables of L 2 ) has its contribution to P(x\, . . . , x n ) 
bounded by 2 C (■^■) D a l ~ 1 . Hence the total contribution of P 2 is bounded in absolute 
value by 2 n+c+Dl ° s( - 1 P'" 1 a' -1 which is equal to a l ~^ . If Pi(x±, . . . ,x n ) is zero, this also 
bounds the absolute value of P(x\, . . . , x n ). 



Assume now that Pi(xi, . . . ,x n ) is not zero, since x L2 is equal to a 1 , the absolute 
value of x L2 P\{x\, . . . , x n ) has a 1 for lower bound. By the triangle inequality 

\P(xi, . . . ,x n )\ > ||a; L2 Pi(a;i,...,x„)| - |P 2 (xi, . . . ,x„)|| 
\P(xi, . . . , x n )\ > a 1 — a l ~ 2 > a l ~ 2 

We can then decide if P\{x\, . . . ,x n ) is zero by comparison of P{x\, . . . , x n ) to 
a'~2 . Remark that Pi(xi, . . . ,x n ) may be zero even if Pi is not zero. Nonetheless P\ 
only depends on variables which are in L\ and are thus randomly taken in [— ]. By 
Lemma [U the probability that the polynomial P\ is not zero although P\(x\, . . . ,x n ) 
has value zero is bounded by e. We have then designed an algorithm which decides 
with probability 1 — e if P has a monomial whose support contains Li and is contained 
in L\ U L2. 

Remark that this implementation of notjzeroJimproved needs only one oracle call 
but requires big evaluation points and to know C in advance. To implement coefficient 
just do the same oracle call and an integer division of the value by a 1 to get the 
coefficient. 

Proof of Theorem [2] : 

The procedure not_zero -improved does one interpolation on a degree I polynomial 
where / is bounded by D. We can bound the value of the polynomial on points of 
value less than D by 2™2 C D D , where C is a bound on the size of the coefficients of 
the polynomial. Since we have assumed that C — 0(n) and that D < n because the 
polynomial is multilinear, the logarithm of the values taken by the polynomial for the 
interpolation is bounded by nlog(n). The univariate interpolation then needs a time 
0(D 2 n\og(n)) and D oracle calls on points of size log(D). 

The procedure not-zero-improved is called in Algorithm [3] with an error parameter 
-4^, it therefore repeats the previously described interpolation 0(n + log(i)) times. 
Each call to not_zero-improved needs a time 0(D 2 n log(n))(n + log(i)) and D(n + 
log(i)) oracle calls. 

Between the visit of two leaves, we call the procedure not_zero -improved at most 
n times and once the procedure coefficient which has a similar complexity. Hence the 
delay is bounded in time by 0(D 2 n 2 log(n))(n + log(i)) and by nD(n + log(i)) oracle 
calls on points of size log(2P). 

Finally since we call the procedures not-zero -improved and coefficient less than nt 
times during the algorithm, the error is bounded by nt^^ < e. □ 

We describe here the way to a improve the error bound for IncPP algorithms, but 
it would work equally well on DelayPP ones. Note that in both cases we need an 
exponential space and there is a slight overhead. 

Proposition 4. // a problem A is in IncPP then there is a polynomial Q and a 
machine M which for all e computes the solution of A with probability 1 — e and satisfies 
for all x, T(x, i + 1) - T(x, i) < Q(\x\,i) log(±). 

Proof. Since A is in IncPP, there is a machine M which computes the solution of 
A with probability | and a delay bounded by Q(\x\,i). Since A(x,y) may be tested 
in polynomial time, we can assume that every output of M is a correct solution, by 
checking A(x,y) before outputting y and stopping if not A(x,y). We now simulate k 



runs in parallel of the machine M on input x. Each time we should output a solution, 
we add it to a set of solution (with no repetition). Assume we have already outputted 
i solutions, we let the k runs be simulated for another Q(\x\,i) steps each before 
outputting a new solution of the set of found solutions and stop if it is empty. This 
algorithm clearly works in incremental polynomial time and if one of the run finds all 
solutions, it also finds all solutions. Then the probability of finding all solutions is more 
than 1 - § . If we set k = ^7§j , we have a probability of 1 — e, which achieves the 
proof. 

Proof of Proposition [2] : 

Assume AnotherSolution^ is computable in probabilistic polynomial time, we 
want to enumerate the solution of the enumeration problem A on the input x. We know 
a bound on the number of solutions of A, that we call B. We assume that the algorithm 
which decides AnotherSolution^ has a probability of error of j^. That is achievable 
by repeating at most log(-B) times the original algorithm, therefore the running time 
is still polynomial. We apply this algorithm to x and the empty set, we add the found 
solution to the set of solutions and we go on like this until we have found all solutions. 
The delay between the i th and the i + 1 th solution is bounded by the execution of the 
algorithm AnotherSolution which is polynomial in \x\ and i the size of the set of 
already found solutions. Moreover the probability of error is bounded by ^ = B x j^. 
This proves that A is in IncPP. 

Conversely if A g IncPP, on an instance (x,S) of AnotherSolution^ we want 
to find a solution which is not in S. We enumerate | *ST | — t- 1 solutions by the IncPP 
algorithm, in time polynomial in \S\ and \x\. If one of these solutions is not in S, it is 
the output of the algorithm. If S is the set of all solutions, the enumeration will end in 
time polynomial in S and x, which allow us to output the value meaning there is no 
other solutions. □ 



